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Abstract 
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In order to overcome the shortcoming of nonconservation of nuclear mo- 
mentum existing in the original Q 2 -rescaling model(OQ 2 RM) and avoid 
■ introducing nuclear shadowing factor, we proposed a double Q 2 -rescaling 

model(D(5 2 RM) for the parton distributions of the bound nucleon. Using 
the experimental data of lepton-nucleus deep inelastic scattering(DIS) and 
the condition of the nuclear momentum conservation, the (5 2 -rescaling pa- 
rameters of various partons for Sn, Fe, Ca and C nuclei are determined. 
The rescaling parameters of valence quark distributions are larger than 
unity and graduauly increases with atomic number A, on the contrary, 
the rescaling parameters of sea quark distributions and gluon distribu- 
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tions are smaller than unity, and slowly decreases with A. By using this 
model, the experimental data of the DIS process, the nuclear Drell-Yan 
process and J/ip photoproduction process are consistently and quite sat- 
isfactorily explained. 

Key words — parton distribution functions, Nuclear effect, Double Q 2 - 
rescaling model. 
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1. Introduction 

In 1982, the European Muon Collaboration at CERN discovered in the 
muon-nucleus deep inelastic scattering(DIS) that the structure functions F 2 of 
bound nucleon is quite different from that of free nucleon within the region 
0.1 < x < 0.7M, this nuclear effect was thus called the EMC effect. Since the 
discovery of the EMC effect, physicists have proposed several models and made 
out more or less satisfactory explanation for it ' 2 '. But, for the shadowing effect 
emerged within very small x region in DIS process ^ and the A-dependence of 
the differential cross section of h-A Drell-Yan process only the constituent 
quark model, Q 2 - rescaling model ^ and the extended x-rescaling model $\ after 
introducing the nuclear shadowing factor, can make out satisfactory explanation. 
In addition, only the three models mentioned above have made out quantitative 
description for J/ip photoproduction on the nuclear target. 
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In the Q 2 -rescaling model with shadowing factor, except for the pa- 
rameter relative to shadowing,for the distribution functions of valence quark, sea 
quark and gluon within the bound nucleon, a common Q 2 -rescaling parameters 
£ (A) was employed. Although the model can explain several kinds experimental 
processes mentioned above, but as pointed by LI Guanglie et. aU 6 ' it leads to 
non-conservation of nuclear momentum, that is, the momentum of a nucleus is 
no longer equal to the sum of momenta of valence quarks, sea quarks and gluons 
including in the nucleus. They thought that for the valence quark distribution 
and sea quark distribution, two different Q 2 - resca lrng parameters should be em- 
ployed. 

As a phenomenological approach to investigate the parton distribu- 
tion functions of a bound nucleon, here we presented a model, in which three 
<5 2 -rescaling parameters were employed for the distributions of valence quark, 
sea quark and gluon respectively, and the artificial shadowing factor was not 
used. Due to restriction of nuclear momentum conservation, only two among the 
three parameters are independent, so we call this model the double Q 2 -rescaling 
model(DQ 2 RM). 

From large amounts of available experimental data, it can be seen that 
for the both of free nucleon and bound nucleon, the x distribution and its evo- 
lution with Q 2 of valence quarks are quite different from those of sea quarks 
and gluons, The distribution of valence quarks is mainly restricted within the 
region of x > 0.1, the " hardness" of the distribution allows one to imagine 
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that the valence quarks are confined within a bag with a distinct boundary and 
the behavior of its evolution is described by GLAP equation The EMC ef- 
fect mentioned above mainly reflects the change of valence quark distribution of 
bound nucleon as compared with free nucleon. For the distributions of sea quark 
and gluon, however, the case is much complex. The distribution of sea quark and 
gluon mainly restricted within the region of x < 0.1, and quickly increases with x 
decreasing. The case is more obvious for the gluon distribution. Therefore, they 
play dominant role in the region of very small x' 8 '. As well known, the longitu- 
dinal size of parton with momentum fraction x is Az = p is the momentum 
of nucleon in a certain frame of reference with high speed, then usual picture 
that partons are confined in a colour bag becomes indistinct. In relation with 
the case, soft partons with very small x value from adjacent nucleons lying in 
the same longitudinal line within nucleus may penetrate and overlap each other. 
The smaller x value ,and the larger the atomic number A, the more outstanding 
the phenomena. Then, recombination or fusion between two or several partons 
from different nucleons could occur. The partons with very small x recombine 
and form into a parton with slightly larger x, and this process naturally leads 
to the shadowing in very small x region and antishadowing in slightly larger x 
region, thereby, after taking into account the recombination, the artificial shad- 
owing factor is not needed. 

Based on above mentioned understanding, we present DQ 2 RM. Using 
this model, the influence of nuclear environment on the parton distributions and 
the experimental data on 1-A DIS process, h-A Drell-Yan process and 7*-A J/ip 
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production can be consistently and satisfactorily explained. In contrast with the 
rescaling parameter £v(A) of valence quark distribution is larger than unity and 
slowly increases with A increasing , a characteristic feature of our model is that 
both the rescaling parameters of sea quark distribution and £g(A) of gluon 

distribution are smaller than unity and slowly decrease with A increasing. This 
feature just right meets the needs of the Q 2 evolution picture of parton distribu- 
tions in small x region of bound nucleon and nuclear momentum conservation. 
We shall present experimental support for this feature in the section 3. 

2. The double Q 2 -rescaling model 
and relative formulas 

For a nucleus A, our model assumes three different Q 2 -rescaling parame- 
ters for the distributions of valence quarks, sea quarks and gluons in the bound 
nucleon respectively: 

q^(x,Q 2 ) = q^(x,^Q 2 ), 

4(*,Q 2 ) = fe^Q 2 ), 
^Q 2 ) = q^(x^iQ 2 ), 

G A (x,Q 2 )=G N (x,^Q 2 ). (1) 

where i denotes the flavor of quark. That is to say, the distributions of various 
partons of bound nucleon may be obtained by Q 2 rescaling for the distributions 
of corresponding partons of free nucleon. For all of nuclei and any Q 2 values, the 
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momentum conservation condition 

1 

/ xdx{Y\qUx, Q 2 ) + <&(x, Q 2 ) + qi(x, Q 2 )] + G A (x, Q 2 )} = 1. (2) 
o * 

should be respected. Because of the restriction of formula(2), only two among 
three Q 2 -rescaling parameters are independent, so this model is called the double 
Q 2 -rescaling model. 

In the DQ 2 RM, the structure function and gluon distribution function 
of bound nucleon respectively is 

F 2 A (x, Q 2 ) = £ e 2 x[q^(x, £ A Q 2 ) + q»(x, £ A Q 2 ) + q^x, ^Q 2 )] (3) 

i 

and 

xG A (x,Q 2 )=xG N (x,£ A Q 2 ), (4) 

F A (x, Q 2 ) denotes the average nucleon structure function of a ideal nucleus with 
the equal number of protons and neutrons(N=Z=|A). 

The nuclear effect of the nucleon structure function and gluon distri- 
bution may be represented by following ratios: 

A (x, Q 2 
F 2 D (x,Q 2 ) 



H A/D (r O 2 ) - (r\ 
G D (x, Q 2 ) ' {b) 

usually, these ratios calculated by a certain theoretical model are used to com- 
pare with the experimental data. 



6 



For the Drell-Yan process in p-A collision, the nuclear effect is repre- 
sented by the ratio 

r d 2 aP- A (x,x t ,Q 2 ) dx 

rpA/N ( -o2\ _ I dxdx t /y\ 

I r d 2 <jP- N (x,x t ,Q 2 ) ^ V) 
J dxdxt 

where 

d2aP ~7X t,Q2) = km* E eU[qUx ' Q2) + qUx ' Q2)] ^ Q2) 

+ f Sl (x,Q 2 )[q^(x t ,Q 2 )+qi(x t ,Q 2 )]}, (8) 

— dxix~ t — = sm^^ eM ^ Q)+qs ^ Q )]q ^ Q ] 

+ qUx,Q 2 M l (x u Q 2 )+q§ t (xt,Q 2 )]}, (9) 
where M 2 denotes the invarint mass of lepton-pair(///i) produced in the process. 



Three <3 2 -rescaling parameters £y, £s and£c are determined by the fol- 
lowing method: firstly, we choose suitable £y and£s to explain the nuclear effect 
in 1-A DIS process; secondly, we describe p-A Drell-Yan process with the same 
values of £y and to test them, and finally determine £g using nuclear momen- 
tum conservation and explain the nuclear effect of J / ip photoproduction process 
at nuclear target. 



For the parton distribution function of free nucleon, we adopt the par- 
ton momentum distributions given by M.Gluck, E.Reya and A.Vogt(GRV)^, 
because this formalism allows to be evoluted inversely to Q 2 ~ l.OGeV 2 . This 
is quite important point for getting the distribution functions of sea quark and 
gluon of bound nucleon by formula(l) at £5 < 1 and £q < I. 
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At present, the statistics and the precision of experimental data on 
1-A DIS process, h-A Drell-Yan process and 7* — A J/ip production process is 
insufficient to test very precise theoretical calculation, so that in our calculation, 
only the leading order approximation on the formulas of hard subprocesses and 
parton distribution functions are considered. 

The Q 2 -rescaling parameters £y and £q on nuclei C 12 , Ca 40 , Fe 56 
and Sn 119 obtained by above method from the data of references[3, 4, 10] are 
listed in the table 1. It is worth to note that, all of are larger than unity 
and increase with A increasing, but all of £g and £q are smaller than unity and 
decrease with A increasing. In the next section , we shall present a experimental 
evidence supporting this result. 



Table. 1 











c 


1.3 


0.7 


0.86 


Ca 


1.35 


0.67 


0.81 


Fe 


1.41 


0.62 


0.76 


Sn 


1.57 


0.45 


0.58 



3. Results and discussion 

Using the values of Q 2 -rescaling parameters listed in the table 1 and formu- 
las (1), (3) and (4), in which LO parton momentum distributions of free nucleon 
given by GRV are adopted' 9 ], we obtained various parton distribution functions 
and average structure function F 2 of bound nucleon in nuclei C 12 , Ca 40 , Fe 56 
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and Sn 119 . Furthermore, using formulas (5)- (9), the theoretical predictions on 
R A l D (x,Q 2 ) ,T A ^ N {x t ,Q 2 ) and Rq /d {x,Q 2 ) by the double Q 2 -rescaling model 
are calculated. 

Theortical curves of R A / D (x,Q 2 ) on 1-A DIS processes for nuclei C 12 , 
Ca 40 , Fe 56 , Sn 119 within the Q 2 range given in reference[3] are separately plotted 
in Fig.l, where corresponding experimental data[3] are shown. It can be seen 
that by our model, only taking different parameters for valence quark distribu- 
tion and sea quark distribution, without introducing the the shadowing factor, 
the experimental data of nuclear effect in 1-A DIS process can be explained quite 
well in the almost whole x region. 

Theoretical curves of T A / D (x t) Q 2 ) on p-A Drell-Yan process for nuclei 
C 12 , Ca 40 and Fe 56 within region of 0.025 < x t < 0.30 and 4 < M fl < 9GeV 
and at Eqm = 4:0GeV are separately plotted in Fig.2a-c, where corresponding 
experimental data' 4 ! are shown. We can see that our results are in qualitative 
agreement with the data, and especially in the region of Xt < 0.15, the agreement 
is quite satisfactory. In addition, within the region 0.15 < x t < 0.30, the agree- 
ment case for heavier nucleus is clearly better than for lighter one. Furthermore 
comparing our results which the prediction given by original Q 2 -rescaling model, 
it can be seen that the DQ 2 RM overcame the shortcoming that the tendency of 
the curve given by original Q 2 -rescaling model is wholly different from the data. 
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Theoretical curve of Rq (x, Q 2 ) (the ratio of gluon distribution func- 
tion of nucleus Sn to one of nucleus C) is plotted in Fig. 3, where the experimental 
data' 10 ' of the ratio measured in 7*+A — > + X process are shown, the result 
calculated by original Q 2 -rescaling model is also plotted in the Fig. 3. Although 
there are quite large errors in the experimental data, it can still be seen that the 
prediction given by our model is clearly better than OQ 2 RM. 

Curves of R$ D {x,Q 2 ) from DQ 2 RM within region 0.1 < x < 0.2 for 
nuclei C,Ca,Fe and Sn are plotted in Fig.4. It can be seen that slope of these 
curves slowly increasing with A increasing. This case reflects that the influence 
of gluon recombination is enhanced with A increasing. 

In the following paragraphs, we shall discuss the problem of the rescal- 
ing parameters £y > l,£s < 1 and £c < 1- 

From above comparision of the calculating results obtained by our 
model with the experimental data we can see that our model has certain ratio- 
nality, and specific values of £y,£s and £g listed in the table 1 are supported by 
the relevant experiments . 

Firstly, in our model, £ v > 1- This is the same as the OQ 2 RM. The 
latter model can well explain the nuclear effect of the DIS process in the region 
0.2 < x < 0.7, where the momentum distribution of nucleus is mainly determined 
by valence quark, so the value of £y in our model should roughly be consistent 
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with the latter. In physics, £y > 1 can be explained by the swelling of nucleon. 
When the nucleon is swollen, according to uncertainty principle, the momentum 
distribution of valence quarks moves to small x region, i.e. the momentum dis- 
tribution of valence quarks was softened. 

Secondly, because of the following reasons , we assume that < 1. 
The first reason is that the OQ 2 RM can not explain the experimental data of 
nuclear effect of the 1-A DIS process and nuclear Drell-Yan process in small x 
region, where the nuclear effect is mainly determined by the change of sea quark 
distribution, so we need to reconsider the Q 2 -rescaling parameter of sea quark. 
The second reason is that some authors of this paper calculated the values of 
the nuclear effect function Rg(x,Q 2 ) of sea quark distribution by using the ex- 
perimental data of 1-A DIS process and nuclear Drell-Yan process^ 11 ' and found 
that in region x < 0.2 R$(x, Q 2 ) increases from the value smaller than unity to 
the value slightly larger than unity with x increasing, and the larger the number 
A, the smaller the value of Rg(x,Q 2 ). If we use Q 2 -rescaling model to describe 
this change of sea quark distribution influenced by nuclear environment, it must 
demands £s < 1, and the larger the A, the smaller the 

Thirdly, we shall present another experimental surport on £q < 1. It is 
well known that the gluon distribution of nucleon is the most important in the 
small x region, and in the region of x < 10~ 2 , gluon distribution almost saturates 
the parton distributions of nucleon. In the case, structure function F 2 (x, Q 2 ) is 
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fully determined by the gluon distribution, and it can be obtained from PQCD^ 12 ! 

dk 2 n .x.a s {k 2 )d[zG{z,k 2 )} 



1 Q 2 



F ^ = ^j dz j-P qA X -) ^ d{lnk2) 



(10) 

where G — > ^ splitting function P g . G (x) = \[x 2 + (1 — a;) 2 ]. Due to the small the 
x value, the more important the gluon distribution, therefore equation (10) can 
be approximately expressed as 

F 2 (x,Q 2 ) ~ ±J2 e * a s(Q 2 )G(x,Q 2 )Ax. (11) 
8tt . 

Recently, NA37/NMC collaboration at CERN measured the ratio F$ n /Ff in the 
region of 0.0125 < x < 0.55 and 1 < Q 2 < 10 2 Ge\/ 2 [ 13 l their result shows that in 
the region of x < 3 x 10~ 2 , Ft f n (x, Q 2 )/ F§ (x, Q 2 ) is clearly less than unity, and 
this result means G Sn (x,Q 2 ) < G c (x ) Q 2 ) in the region. If we use <5 2 -rescaling 
approach to get G Sn (x, Q 2 ) from G c (x, Q 2 ) and assume 

G Sn (x,Q 2 )=G c \x,^ c Q 2 ), (12) 

where £ Sn / c should be smaller than unity ,that is to say, the gluon distribution 
function of a heavier nucleus may be obtained from one of a slighter nucleus 
by Q 2 rescaling approach and using Q 2 rescaling parameter smaller than unity. 
Therefore, we may conjecture that the gluon distribution functions of all nuclei 
may be obtained from one of free nucleon or deuterium by this approach and 
using £q smaller than unity, moreover, the values of £^ slowly decreases with 
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atomic number A increasing. 

In summary, we presented DQ 2 RM, in which different rescaling param- 
eters £y,£s and £g are assumed, and the shadowing factor is not introduced. By 
this model, the experimental data of the 1-A DIS process, the nuclear Drell-Yan 
process and 7* — A J/ip production process were consistently and quite satisfac- 
torily explained. But, this model is a initial and phenomenological model, on 
the one hand, it needs to be further improved as more precise experimental data 
presented, on the other hand, its physical foundation still needs to be explored 
more elaboratively. 

The authors are grateful to Dr. Cong- Feng QIAO for useful help in com- 
puting work. 
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Figure Captions 
Fig.l The ratio R A ' D (x,Q 2 ) for nuclei C 12 , Ca 40 , Fe 56 , Sn 119 versus 
the momentum fraction x within the Q 2 range given in reference [3] , where cor- 
responding experimental data [3] are shown. The solid lines are the results of 
our model, the dashed lines are the results of the old Q 2 — rescaling model(In the 
OQRM, the Q 2 -rescaling parameters £ for nuclei C 12 , Ca A0 , Fe 56 , Sn 119 are 1.6, 
1.86, 2.02, 2.24, respectively). 

Fig.2 The ratio T A / D (x t ,Q 2 ) for nuclei C 12 , Ca m and Fe 56 versus x 
within region of 0.025 < x t < 0.30 and 4 < M a < 9GeV and at E CM = 40GeV, 
where corresponding experimental data' 4 ' are shown. The meaning of lines is the 
same as that in Fig.l. 

Fig.3 The ratio R S G n,C \x, Q 2 ) versus x, where the experimental data' 10 ! 
are shown. The meaning of lines is the same as that in Fig.l. 

Fig.4 The ratio Rg /D (x,Q 2 ) from DQ 2 RM versus x within region 
0.1 < x < 0.2 for nuclei C 12 , Ca m and Fe 56 . 
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